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I . Numerical  Methods  for  Control  and  Identification 
of  Delay  Systems  (H.T.  Banks) 


1 . Spline-based  techniques 

In  the  past  year  Banks  and  Kappel  (Universitat  Graz, 

Austria)  have  developed  ideas  involving  the  use  of  <=p!  ine-based 
approximations  for  computation  of  solutions  of  systems  of  equa- 
tions with  delays.  Through  their  theoretical  and  numerical 
investigations,  they  have  strong  evidence  to  support  the  use  of 
these  approximations  as  an  attractive  alternative  to  the  "averaging" 
type  approximations  (see  [1])  in  optimal  control  and  para- 
meter identification  problems.  The  basic  conceptual  framework 
for  the  spline  methods  along  with  a summary  of  numerical  results 
for  integration  of  linear  delay  systems  can  be  found  in  a just 
completed  manuscript  [2].  In  these  efforts,  numerical  experiments 
using  both  piecewise  linear  and  piecewise  cubic  spline  bases  were 
carried  out  with  excellent  results.  Pursuing  these  ideas.  Banks 
developed  software  packages  to  use  the  spline  methods  in  optimal 
control  problems.  Numerical  experiments  were  made  with  test 
examples  and  again  it  was  found  that  the  spline-based  methods 
appear  to  offer  significant  advantage  over  the  "averaging"  methods 
for  control  problems.  Details  of  some  of  these  numerical  findings 
can  be  found  in  [3] . Efforts  (joint  with  J.  Burns  and  E.  Cliff, 
an  aerospace  engineer  at  VPI)  have  now  begun  to  use  these  ideas 
in  the  context  of  parameter  estimation  schemes.  The  early  results, 
both  theoretical  and  numerical,  are  extremely  promising. 
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. Approximation  of  nonlinear  systems 


In  [4]  Banks  developed  a theory  based  on  the  averaging 
approximations  to  treat  nonlinear  delay  system  optimal  control  problems, 
Numerical  results  to  demonstrate  the  feasibility  of  these  methods 
were  also  given.  Banks  has  now  made  some  progress  in  extending 
these  ideas  to  allow  for  more  general  nonlinear  systems  and  also 
to  allow  for  spline  and  other  type  of  approximation  schemes  as  well 
as  the  "averaging"  schemes.  The  theory,  when  completed,  will  hope- 
fully also  be  applicable  to  nonlinear  parameter  estimation  problems. 
Efforts  on  both  theoretical  and  numerical  aspects  of  the  work  is 
continuing . 


3 . Difference  equation  type  techniques 

D.  Reber,  under  the  direction  of  Banks,  has  completed  his 
Ph.D.  thesis  [5]  on  modifications  (of  the  Banks-Burns  theoretical 
framework  [1])  which  result  in  difference  equation  approximations 
(as  opposed  to  the  finite  systems  of  ordinary  differential  equa- 
tions in  the  Banks-Burns  approach) . Both  theoretical  and 
computational  efforts  to  study  general  linear  nonautonomous 
(time-varying)  system  control  problems  (via  the  "averaging"  type 
approximation  scheme)  were  made.  The  results  (which  are  summarized 
in  the  recent  manuscript  [6])  indicate  that  while  the  methods 
developed  by  Reber  offer  a viable  alternative,  one  cannot  expect 
to  make  significant  gains  here  unless  one  employs  something 
better  than  a simple  first  order  differencing  scheme  on  derivatives 
in  the  problems.  Banks,  and  another  graduate  student,  Rosen,  are 
now  trying  to  develop  ideas  to  handle  higher-order  differencing 
approximations.  When  combined  with  the  higher-order  spline  based 
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schemes,  these  ideas  should  yield  methods  considerably  more 
efficient  than  those  develpped  to  date. 


I I . Optimal  Control  of  Diffusion-Reaction  Systems  (H.T.  Banks) 
Banks,  J.P.  Kernevez  and  M.  Duban  (an  assistant  to  Kernevez 
at  Universite  de  Technologie  de  Compiegne)  are  continuing  their 
efforts  on  methods  for  optimization  of  diffusion-reaction  problems 


ds  d s ^ a ✓ 
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through  boundary  controls.  (Here  v is  a nonlinear  reaction 
velocity  approximation  determined  by  the  specific  reactions 
considered.)  Some  progress  in  the  theoretical  aspects  of  these 
investigations  has  been  made  and  is  partially  summarized  in  [7] . 
Additional  theoretical  results  have  been  obtained  which  reveal  to 
some  extent  how  "singularities"  arise  in  the  standard  numerical 
procedures  (i.e.  conjugate-gradient  methods)  when  applied  to  these 
problems  even  though  the  problem  itself  may  be  theoretically  well- 
posed.  Current  efforts  are  centered  on  comparison  of  several 
alternative  procedures  for  efficiently  finding  numerical  solutions, 


III.  Qualitative  Theory  of  Delay  Systems  (J.K.  Hale) 

1 , Effects  of  variations  in  the  delays 

In  many  problems  in  the  applications,  one  encounters 


differential  difference  equations 
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(1)  x(t)  = f (x(t) ,x(t-r1) , . . . ,x(t-rn)) 
and  difference  equations 

(2)  y(t)  = g(y(t-rj) . fy(t-rn)) 

where  0 < r^  < r^  < ...  < r . The  vector  parameter  r = (rj,...,rn) 
is  a physical  parameter  which  is  known  only  to  belong  to  some 
interval.  It  is,  therefore,  important  to  know  how  properties  of 
the  solutions  depend  upon  r. 

Even  for  the  difference  equation  (2),  this  problem  must  be 
investigated  for  initial  data  from  a function  space  & of  func- 
tions on  [-b,0]  where  b > r . Hale  and  some  of  his  students 
have  been  investigating  this  question  for  some  time.  The  basic 
difficulty  centers  around  the  fact  that  the  solution  x(<J>,r) 
through  an  initial  function  can  never  be  smooth  in  r regard- 
less of  the  function  space  The  overall  objective  is  to 

determine  a qualitative  theory  even  though  these  smoothness 
properties  in  r do  not  hold. 

Specific  results  obtained  in  the  past  year  are  the  following. 
For  a linear  difference  equation 

(3)  y (t)  = l A . y (t - r . ) 

j=l  J J 

where  each  Aj  is  an  N * N constant  matrix,  the  asymptotic  be- 
havior of  the  solutions  is  given  by  the  set 

-Xr 

Z(r)  = (Re  X:  det[I  - E A ^ e J]  = 0}.  Hale  and  Avellar  have  proved 
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the  following: 

Theorem.  v There  exist  continuous  functions  a(r),p(r)  such  that 

cl  Z(r)  c [P (r) ,o(r) ] 

and  this  is  the  minimal  interval  containing  cl  Zfr)  if  the 
components  of  r are  rationally  independent. 

A method  is  also  given  for  computing  P(r),o(r).  If  Tg(t,r) 
is  the  semigroup  generated  by  (3)  on  CQ  = {^  e C ( [ -b , 0]  dRn) : 

^ (0)  = E Aj'P  (-tj  ) } defined  by  (TQ (t , r)^)  (0 ) = yO  ,r)  (t+6)  , 

-r  < 0 < 0,  and  y , r ) is  the  solution  of  (3)  through  ^ , then 
they  also  prove  that  Y(Tg(t,r))  = re(TQ(t,r))  < exp(o(r)t)  with 
equality  holding  for  the  components  of  r rationally  independent. 
Here  Y is  the  spectral  radius  and  rg  is  the  radius  of  the 
essential  spectrum. 

The  implications  of  this  result  for  stability  are  clear.  Hale 
and  Avellar  also  show  how  this  result  can  be  used  to  obtain  con- 
tinuous dependence  on  r of  solutions  of  the  neutral  equation 

[x(t)  - E AjX(t-Tj ) ] = BQx (t)  + E BjX(t-r^ ) . 

These  results  have  appeared  in  the  thesis  of  Avellar  and  are  being 
prepared  for  publication. 

In  an  attempt  to  discover  the  essential  elements  of  a 


qualitative  theory  for  (1)  as  a function  of  r,  Hale  has  discovered 
how  to  prove  that  a smooth  Hopf  bifurcation  with  respect  to  r 
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occurs  for  (1)  even  though  the  solutions  are  not  differentiable 
in  r ([8]).  This  has  important  implications  to  the  existence 
of  periodic  solutions  for  equations  with  several  delays.  The 
methods  should  be  applicable  to  other  problems. 

2 . Dissipative  systems 

Cooperman  [9],  a student  of  Hale,  has  completed  his 
dissertation  on  a-contractions  and  dissipative  processes.  The 
results  have  implications  to  the  asymptotic  behavior  of  the  solu- 
tions for  retarded  and  neutral  functional  differential  equations. 
In  addition,  he  has  shown  an  intimate  connection  between 
dissipative  systems  and  Liapunov  functions. 

IV.  Nonlinear  Oscillations  and  Bifurcation  Theory  (J.K.  Hale) 

1 . Multiple  parameter  bifurcation  problems 

Hale  and  Taboas  [10,  11]  have  completed  their  work  on 
harmonic  and  subharmonic  bifurcation  in  the  equation 

x + G (x)  = Ax  + pg(t) 

where  k,P  are  small  independent  parameters,  g(t+l)  = g(t)  and 
the  nonlinear  equation  x + G(x)  = 0 has  a nonconstant  periodic 
solution  of  period  k,  an  integer.  The  methods  have  new  applica- 
tions to  homoclinic  orbits,  but  the  results  are  not  complete  at 
this  time. 

2 . Reaction-diffusion  equations 

Alikakos  [12,  13]  a student  of  Hale  and  Mallet-Paret , has 
completed  his  dissertation  on  asymptotic  behavior  in  parabolic 


equations.  He  obtained  results  on  a scalar  parabolic  equation 
with  nonlinear  boundary  conditions,  generalizing  results  of  Aronson, 
Peletier  and  Ball.  He  also  obtained  results  on  systems  of  equations 
with  food  pyramid  conditions,  generalizing  results  of  Chow  and 
Williams.  In  addition  to  showing  that  the  methods  of  dynamical 
systems  are  more  powerful  than  maximum  principle  arguments,  he 
was  able  to  obtain  uniform  bounds  on  the  solutions  rather  than 
bounds  with  1 < p < 00 . 

V.  Control  Systems  (J.P.  LaSalle) 

1 . Stable  feedback  generators 

An  initial  effort  on  this  problem  by  LaSalle  and  Palmer 
has  been  discontinued.  Some  progress  was  made  on  attempting  to 
extend  the  concept  of  solutions  of  differential  equations  with 
discontinuous  right-hand  sides.  Earlier  work  of  Fillipov,  Hermes, 
Hajek,  and  others  was  extended.  However,  recent  work  by  Aizerman 
and  one  of  his  colleagues  at  the  Institute  of  Control  Sciences, 
Moscow,  in  this  same  direction  indicates  it  does  not  provide  the 
necessary  mathematical  framework  for  the  problem  of  investigating 
stable  feedback  generators. 

2 . Learning,  identification  and  adaptation 

The  concept  of  eventual  stability  for  nonautonomous  ordinary 
differential  equations  was  introduced  by  LaSalle  around  1958  and 
studied  by  R.  Rath  in  his  Notre  Dame  Ph.D.  dissertation,  1962.  An 
application  was  made  in  1963  (R.  Rath  and  J.P.  LaSalle,  "Eventual 
Stability",  Proceedings  of  Second  IFAC  Congress,  Basel,  1963)  to 
prove  the  convergence  of  a scheme  for  adaptive  control.  It  is  clear 
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today  that  the  concept  of  eventual  stability  is,  and  should  be, 
related  to  the  limiting  equations  of  nonautonomous  systems.  Last 
fall,  LaSalle  and  Zvi  Artstein,  Weizmann  Institute  of  Science, 

Israel,  began  a joint  investigation  of  the  relationship  of  ques- 
tions of  stability,  adaptation,  identification,  and  learning  to 
the  limiting  equations  of  control  systems.  This  research  is  into 
its  initial  stages,  and  it  is  hoped  that  Artstein  can  visit 
Brown  next  spring  to  continue  these  efforts. 

3 . Stability 

a . functional  difference  equations 

Under  the  direction  of  LaSalle,  Palmer  has  completed 
his  thesis  [14]  on  an  extension  of  Liapunov's  direct  method  for 
the  study  of  the  stability  of  nonautonomous  functional  (delay- 
differential)  equations. 

b . stability  of  nonautonomous  difference  equations 

M.  Latina,  R.I.  Jr.  College,  under  the  direction  of 
LaSalle  has  completed  an  extension  of  the  invariance  principle 
and  Liapunov's  direct  method  for  nonautonomous  discrete  systems 
(time-varying  difference  equations).  LaSalle  gave,  in  outline, 
the  development  of  this  theory  in  his  lectures  at  an  NSF-CBMS 
Regional  Conference  on  Applied  Mathematics,  Mississippi  State, 

August  1975.  Improvements  in  the  theory  have  been  made  that 
extend  the  range  of  applications. 

4 . Stability  and  control  of  discrete  processes 

LaSalle  spent  most  of  his  time  this  summer  working  on  a 
two-volume  exposition  on  the  control  and  stability  of  discrete-time 
processes.  A final  draft  of  Volume  I on  linear  systems  is 
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almost  completed.  Volume  II,  which  will  be  directed  towards  the 
more  general  nonlinear  thepry  of  both  stability  and  control,  will 
be  completed  by  the  end  of  this  year.  Much  of  our  knowledge  of 
discrete  systems  appears  within  many  different  contexts  (numerical 
analysis,  engineers'  studies  of  sampled  data  systems,  probability 
and  statistics,  etc.),  and  this  exposition  provides  through  the 
use  of  the  language  and  concepts  of  dynamical  system  both  a unity 
and  extension  of  the  theory  and  is  a good  preparation  for  the  study 
of  cont inuous - t ime  systems  (ordinary  and  functional  differential 
equations , etc . ) . 
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SPUN!:.  APPROXIMATIONS 


1 


FOR 

FUNCTIONAL  DIFFERENTIAL  EQUATIONS 

II.  T.  Ranks  and  F.  Kappel 

Abstract : We  develop  an  approximation  framework  for  linear  hereditary 

systems  which  includes  as  special  cases  approximation  schemes  employing 
splines  of  arbitrary  order.  Numerical  results  for  first  and  third 
order  spline  based  methods  are  presented  and  compared  with  results 
obtained  using  a previously  developed  scheme  based  on  averaging  ideas. 


APPROXIMATION  OF  DELAY  SYSTEMS  WITH  APPLICATIONS 


TO  CONTROL  AND  IDENTIFICATION 


II.  T.  Ranks 


i 


Abstract : We  discuss  approximation  ideas  for  functional  differential 

equations  and  how  these  ideas  can  be  employed  in  optimal  control  and 
parameter  estimation  problems.  Two  specific  schemes  are  described, 
one  based  on  integral  averages  of  the  function  being  approximated, 
the  other  based  on  best  L^  spline  approximations.  An  example 
illustrating  numerical  behavior  of  these  schemes  applied  to  an  optimal 
control  problem  is  presented. 
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OPTIMAL  CONTROL  OF  DIFFUSION-REACTION  SYSTEMS 


by 

II.  T.  Banks,  fl.  C.  Duban  and  J.  P.  Kernevez 

Abstract:  We  discuss  several  formulations  of  optimization 

problems  which  arise  in  a natural  way  in  the  investigation  of 
transport  properties  of  artificial  membranes  and  general  dif- 
fusion-reaction media.  Nonlinear  reaction  velocity  approxima- 
tions dictated  by  reactions  of  interest  to  biochemists  place 
the  problems  in  a class  to  which  one  cannot  apply  the  usual 
computational  techniques  (e.g.  gradient,  conjugate-gradient) 
in  a straightforward  manner.  The  inherent  difficulties,  how 
one  might  circumvent  them,  and  some  of  our  initial  efforts 
towards  development  of  feasible  computational  schemes  are  dis- 


cussed . 


APPROXIMATION  TECHNIQUES  FOP  CONTROL  SYSTEMS  WITH  DELAYS 


by 

H.  T.  Banks  and  J.  A.  Burns 

Abstract 

We  present  a theoretical  framework  for  approximation 
techniques  for  nonlinear  system  optimal  control  problems.  Two 
particular  approximation  schemes  that  may  be  used  in  the  context 
of  this  framework  are  discussed  and  typical  numerical  results  for 
two  examples  to  which  we  have  applied  these  schemes  are  given.  We 
conclude  with  a brief  survey  of  related  investigations. 


HEREDITARY  CONTROL  PROBLEMS:  NUMERICAL  METHODS 


BASED  ON, AVERAGING  APPROXIMATIONS 

by 

H.  T.  Banks  and  J.  A.  Burns 

Abstract:  An  approximation  scheme  involving  approximation  of  linear 

functional  differential  equations  by  systems  of  high  order  ordinary 
differential  equations  is  formulated  and  convergence  is  established 
in  the  context  of  known  results  from  linear  semigroup  theory. 
Applications  to  optimal  control  problems  are  discussed  and  a summary 
of  numerical  results  is  given.  The  paper  is  concluded  with  a brief 
survey  of  previous  literature  on  this  class  of  approximations  for 
systems  with  delays. 


Abstract  of 


a -CONDENSING  MAPS  AND  DISSIPATIVE  SYSTEMS 
by  Gene  David  Cooperman,  Ph.C.,  Brown  University,  June  1978 

Dissipative  systems  are  a generalization  of  differential  equa- 
tions with  Lyapunov  functions.  The  main  idea  is  that  solutions  should 
enter  a given  bounded  set  and  stay  there  for  all  further  time.  It 
was  first  studied  by  N.  Levinson  in  1944.  The  concept  is  intimately 
connected  with  the  existence  of  a maximal,  compact,  invariant  set. 
This  is  important,  since  any  maximal,  compact,  invariant  set  must 
contain  all  equilibrium  solutions,  periodic  solutions,  and  almost 
periodic  solutions.  Its  role  is  analogous  to  the  role  of  the 
maximal,  invariant  set  in  the  LaSalle  Invariance  Principle. 

In  Euclidean  space,  the  various  possible  definitions  of  a 
dissipative  system  are  degenerate.  In  a Banach  space,  this  is  not 
so.  Since  the  natural  setting  for  delay  differential  equations  is 
in  a Banach  space,  it  is  important  to  identify  which  definition  of 
dissipative  systems  should  be  used  in  a particular  application. 
Counterexamples  are  given  to  indicate  where  the  best  possible  result 
has  been  obtained.  Applications  include  Lyapunov  and  stability 
theory,  a continuous  dependence  result  for  invariant  sets,  and  a 
demonstration  of  the  use  of  these  methods  in  facilitating  work  with 
skew  product  flows. 

In  Chapter  II,  a theorem  is  proven  which  allows  us  to  extend  the 
results  on  ct -contractions  to  ot-condensing  maps.  This  is  important, 
since  a -condensing  maps  are  very  general  solution  operators  of  delay 
differential  equations,  which  include  a -contractions.  Additional 
applications  of  the  theorem  include  a generalization  of  a fixed 
point  theorem,  and  of  a continuous  dependence  theorem  for  fixed 
points,  to  the  a -condensing  case. 


NONLINEAR  OSCILLATIONS  IN  EQUATIONS  WITH  DELAYS 


Jack  K.  Male 


Abstract:  These  lectures  arc  concerned  only  with  some  aspects  of 

bifurcation  theory  in  the  local  theory  of  nonlinear  oscillations  in 
equations  with  delays;  that  is,  behavior  of  solutions  near  an 
equilibrium.  In  particular,  we  study  how  the  qualitative  behavior 
of  solutions  change  as  parameters  vary.  A detailed  study  of  the 
local  theory  is  important  in  order  to  know  the  types  of  solutions  to 
expect  in  a global  problem.  Of  course,  there  is  no  reason  to  only 
study  local  theory  near  an  equilibrium.  One  should  study  how  the 
qualitative  behavior  changes  near  any  invariant  set  - for  example, 
behavior  near  a periodic  orbit,  behavior  near  an  orbit  which  connects 
a saddle  point  to  itself,  etc.  More  complicated  behavior  is  expected 
near  these  large  invariant  sets.  One  can  obtain  invariant  torii, 
homoclines  points  which  exhibit  a chaotic  behavior,  etc.  We  restrict 
ourselves  in  these  lectures  to  behavior  near  equilibrium. 

The  simplest  type  of  smooth  bifurcation  is  from  an  equilibrium  to 
a periodic  orbit  - the  so-called  Mopf  bifurcation.  In  Section  2,  we 
discuss  the  Hopf  bifurcation  in  equations  with  finite  delays  per- 
mitting the  bifurcation  parameters  to  be  the  delays  themselves.  At 
first  glance,  such  a result  does  not  seem  possible  because  the  vector 
field  in  the  equation  is  not  differentiable  in  the  parameters.  The 
theorem  does  require  some  new  ideas  and,  for  this  reason,  the  proof 
is  given  in  some  detail.  Several  examples  are  given  in  Section  3. 
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In  Section  4,  similar  results  are  presented  for  equations  with 
infinite  delays.  In  Section  5,  we  give  an  example  in  two  dimensions 
for  which  stable  llopf  bifurcation  occurs  with  decreasing  delay.  In 
Section  6,  we  give  an  introduction  to  some  of  the  methods  available 
for  nonautonomous  equaitons. 


RETARDED  EQUATIONS  WITH  INFINITE  DELAYS 


Jack  K.  Hale 

Abstract;  It  is  the  purpose  of  these  notes  to  describe  the  theory 
of  Hale  and  Kato  for  functional  differential  equations  based  on  a 
space  of  initial  data  which  satisfy  some  very  reasonable  axioms.  We 
also  indicate  some  recent  resvilts  of  Naito  showing  how  extensive 
the  theory  of  linear  systems  can  be  developed  in  an  abstract  setting 
in  particular,  the  characterization  of  the  spectrum  of  the  in- 
finitesimal generator  together  with  the  decomposition  theory  and 
exponential  estimates  of  solutions. 


SOME  RECENT  RESULTS  ON  DISSIPATIVE  PROCESSES 


.lack  K.  Hale 


Abstract : A detailed  study  of  dissipative  processes  has  obvious 

implications  in  the  asymptotic  behavior  of  solutions  of  differential 
equations.  The  properties  of  such  maps  has  been  extensively 
studied  over  the  past  years  with  the  base  space  being  a Banach  space 
rather  than  IRn.  The  results  on  asymptotic  behavior  pertain  to  the 
properties  of  solutions  of  functional  differential  equations  of 
retarded  and  neutral  type  as  well  as  certain  types  of  partial 
differential  equations.  It  is  the  purpose  of  this  lecture  to  present 
some  recent  developments  in  this  theory,  especially  the  ones  due  to 
Coooperman.  A brief  discussion  of  the  history  of  the  subject  may 
be  found  in  Chapter  4 of  Theory  of  Functional  Differential 
Equations  by  J.K.  Hale. 


TOPICS  IN  LOCAL  B' FURCATION  THEORY 


by 

Jack  K.  Hale 

Abstract 

Suppose  A,X,Z  are  Banach  spaces,  M:A  x X -►  Z is  a 
mapping  continuous  together  with  derivatives  up  through  some 
order  r.  A bifurcation  surface  for  the  equation  (1) 

M ( A , x)  =0  is  a surface  in  parameter  space  A for  which 
the  number  of  solutions  x of  (1)  changes  as  A crosses 
this  surface.  Under  certain  generic  hypotheses  on  M,  the 
author  and  his  colleagues  have  shown  that  one  can  systematic- 
ally determine  the  bifurcation  surfaces  by  elementary  scaling 
techniques  and  the  implicit  function  theorem.  This  talk  gives 
a summary  of  these  results  for  the  case  of  bifurcation  near 
an  isolated  solution  or  families  of  solutions  of  the  equation 
M(Ag,x)  = 0.  The  results  have  applications  to  the  buckling 
theory  of  plates  and  shells  under  the  effect  of  external  forces, 
imperfections,  curvature  and  variations  in  shape.  The  results 
on  bifurcation  near  families  has  applications  in  nonlinear 
oscillations  and  the  theory  of  homoclinic  orbits. 


BIFURCATION  NEAR  DEGENERATE  FAMILIES 


by 

Jack  K.  Hale  and  Placido  Taboas 

Abstract : Suppose  X,  A,Z  are  Banach  space,  M:  X x A ->  7,  is  a 

smooth  function  and  the  equation 

M(x,A)  = 0,  (1) 

for  A = 0,  has  a one  parameter  family  of  smooth  solutions  p(s), 

0 < s < 1,  p (0)  = p(l)  with  dp(s)/ds  / 0.  If  Y = (p(s),  0 < s < 1}, 

the  objective  is  to  discuss  the  solutions  of  (1)  near  f for  A in 
a neighborhood  of  zero.  Assuming  A(s)  = 9M (p (s ) , 0) /9s  is 
Fredholm  of  index  zero  and  dim  _y^(A(s))  = 1,  for  0 < s < 1,  the 
authors  have  previously  given  a solution  to  this  problem  under  natural 
hypotheses  on  M.  In  this  paper,  the  case  where  dim  ^flA(s))  = 2 
is  considered.  Applications  are  given  to  the  effect  of  small  damping 
and  small  forcing  on  the  existence  of  periodic  solutions  near  a 

periodic  solution  of  an  autonomous  Hamiltonian  system. 
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BIFURCATION  NEAR  FAMILIES  OF  SOLUTIONS 

-Jack  K.  Hale 

Abstract : Many  investigations  in  bifurcation  theory  are  concerned 

with  the  following  problem.  If  M(0,0)  = 0 and  3M(0,0)/3x  has 
a nontrivial  null  space,  find  all  solutions  of  the  equation 

M (x , X ) = 0 (1.1) 

for  (x,X)  in  a neighborhood  of  (0,0)  r.  x * A. 

If  dim  A = 1;  that  is,  there  is  only  one  parameter  involved 
then  the  existence  of  more  than  one  solution  in  a neighborhood  of 
zero  can  be  proved  by  making  assumptions  only  about  3M(0,0)/3x  and 
3M(0,0)/3x3X  . However,  if  dim  A >2,  then  the  problem  is  much  more 
difficult  and  more  detailed  information  is  needed  about  the  function 
M.  A careful  examination  of  the  existing  literature  for  dim  A > 2 
reveals  that  the  additional  conditions  imposed  on  M imply,  in 
particular,  that  the  solution  x = 0 of  the  equation 

M (x, 0)  = 0 (1.2) 

is  isolated  (see,  for  example,  the  papers  on  catastrophe  theory) . 

These  hypotheses  eliminate  the  possibility  that  Equation  (1.2)  has 
a family  of  solutions  containing  x = 0.  Such  a situation  occurs, 
for  example,  for  M(x,X)  = Ax  + N(x,X),  where  A is  linear  with  a 
nontrivial  null  space  and  N(x,0)  = 0 for  all  x.  There  also  are 
interesting  applications  where  Equation  (1.2)  is  nonlinear  and  there 
exists  a family  of  solutions.  For  example,  Equation  (1.2)  could  be 
an  autonomous  ordinary  differential  equation  with  a nonconstant 
periodic  orbit  of  period  2 n with  the  family  of  solutions  being 
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obtained  by  a phase  shift.  -When  the  differential  equation  in  the 
latter  situation  is  a Hamiltonian  system,  the  parameters  (A^,A2) 
could  correspond  to  a small  damping  term  and  a small  forcing  term  of 
period  211.  To  the  author's  knowledge,  the  first  complete  investiga- 
tions of  special  problems  of  each  of  these  latter  types  are  contained 
in  papers  by  Hale,  Taboas  and  Rodrigues. 

It  is  the  purpose  of  this  paper  to  begin  the  investigation  of 
the  abstract  problem  for  Equation  (1.1),  especially  to  extend  the 
results  in  the  paper  by  Hale  and  Taboas. 
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INTERACTION  OF  DAMPING  AND  FORCING  IN  A SECOND  ORDER  EQUATION 


by 

Jack  K.  Hale  and  Placido  Tciboas 

Synopsis . Suppose  A,y  are  real  parameters,  f is  a scalar  func- 
tion which  is  2Tr-periodic,  xg(x)  > 0 for  x ^ 0 and  consider  the 
equation 

x + g(x)  = -Ax  + yf(t).  (1) 

For  A = y = 0,  every  solution  has  the  form  x(t)  = 4>  (tu  (a ) t + a,  a) 
for  some  constants  a, a and  <p  (6+2n  ,a)  = <|>(0,a).  If  there  is  an 
aQ  such  that  <o  (a^)  = 1 (i.e.,  there  is  a 2Tr-periodic  orbit  T 

in  (x,x)-space)  and  a)'(ag)  / 0,  the  problem  is  to  characterize 
the  number  of  2ir-periodic  solutions  of  Equation  (1)  which  lie  in 
a neighborhood  of  r for  (A,y)  in  a small  neighborhood  of  (0,0). 

A complete  solution  of  this  problem  is  given  under  the  hypothesis  that 
the  function  h(a)  = 

derivative  except  at  a finite  number  of  points  ou  and  h"  (ou)  / 0. 
The  bifurcation  curves  in  (A, y) -space  are  determined  by  the  ct ^ 
and  are  tangent  to  the  straight  lines  A = h(otj)y  at  (A,y)  = (0,0). 
In  general,  the  2ir-periodic  solutions  of  (1)  are  not  continuous  at 
(A»y)  = (0,0).  The  nature  of  this  discontinuity  is  discussed  in 
detail.  it  is  also  shown  that  a necessary  and  sufficient  condition 
for  a 2ir-periodic  solution  x(A,y)  to  be  continuous  at  (A,y)  = (0,0) 
is  that  A/y  constant  as  A -*■  0,  y -*■  0. 
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GENERIC  BIFURCATION  WITH  APPLICATIONS 


J.  K.  Hale 


Abstract : This  paper  is  a set  of  lecture  notes  on  generic 

bifurcation  and  its  applications  with  the  emphasis  on  equations 
involving  more  than  one  independent  parameter.  The  general 
theory  is  discussed  for  problems  which  are  degenerate  to  order 
one  or  two.  Applications  are  given  primarily  to  the  buckling 
of  plates  and  shells  with  the  parameters  representing  external 
forces,  loading,  imperfections,  curvature  and  dimension. 


RESTRICTED  GENERIC  BIFURCATION 


Jack  K.  Hale 


Abstract : In  the  past  few  years,  there  lias  been  considerable  attention 

devoted  to  the  existence  of  bifurcation  for  one-parameter  families  of 
mappings.  Concurrent  with  this  development  has  been  the  extensive 
theory  of  universal  unfolding  of  mappings  or  generic  bifurcation  for 
families  of  mappings  which  depend  on  a sufficiently  large  number  of 
parameters.  The  purpose  of  this  paper  is  to  discuss  methods  for 
determining  the  nature  of  bifurcation  when  the  family  of  mappings  has 
k > 1 parameters  but  k is  generally  smaller  than  the  number  of 
parameters  necessary  to  describe  the  universal  unfolding. 
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NEW  STABILITY  RESULTS  FOR  NONAUTONOMOUS  SYSTEMS 


J.P.  LaSalle 

Abstract 

The  new  invariance  properties  that  have  been  established 
for  nonautonomous  ordinary  differential  equations  greatly  extend 
the  range  and  power  of  Liapunov’s  direct  method  for  the  study  of 
the  stability  of  time-varying  systems.  An  essential  feature  of 
the  method  is  the  establishment  of  a relationship  between  Liapunov 
functions  and  the  location  of  the  positive  limit  sets  of  solutions. 
The  principal  contribution  of  this  paper  is  a theorem  connecting 
Liapunov  functions  and  positive  limit  sets  of  sufficient  generality 
to  close  a gap  in  the  present  theory. 
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STABILITY  THEORY  FOR  DIFFERENCE  EQUATIONS 


by 

J.  P.  LaSalle 

Abstract : This  article  is  designed  to  give  through  the  study  of 

difference  equation  (discrete  dynamical  systems)  a view  of  and  an 
introduction  to  the  general  theory  of  the  stability  of  dynamical 
systems  in  its  most  modern  aspect.  Much  of  what  is  presented 
here  is  known,  although  not  perhaps  as  well  known  as  it  should  be, 
and  there  are  some  things  that  are  new.  One  of  these  has  to  do 
with  a connectedness  property  of  the  positive  limit  sets  of  the 
solutions  of  difference  equations  which  provides  a means  through 
the  use  of  Liapunov  functions  of  establishing  the  existence  of 
equilibrium  points  (fixed  points)  and  oscillations  (periodic  points). 
Another  is  the  generalization  of  the  usual  concept  of  a vector 
Liapunov  function,  and  this  leads  to  a possible  method  of  designing 
control  systems  where  the  measure  of  the  error  or  the  performance 
criterion  is  a vector  rather  than  a scalar.  Applications  of  the 
theory  are  illustrated  by  simple  examples. 

The  article  was  written  for  undergraduate  teachers  of 


mathematics  but  it  should  also  serve  as  a good  introduction  for 
engineers  and  scientists  to  the  latest  results  in  the  theory  of 
the  stability  of  dynamical  systems. 


LIAPUNOV  STABILITY  THEORY  FOR 
NONAUTONOMOUS  FUNCTIONAL  DIFFERENTIAL  EQUATIONS 

John  Wrlliam  Palmer 


A}  >stract 

T t.  i s a standard  technique  in  all  areas  of  ma l.hoina  1 i.c:s  t:o 
reduce  a more  comp  1 i on  hod  problem  I o one  that  has  already  boon 
solved.  Aside  from  providinq  a mi  i ty  to  mathematics,  this  idea  pre- 
sents us  with  new  approaches  !:■'  unsolved  proh l oms . The  area  of 
differential  aquations  lias  man,  interesting  upal  i en(  ions  o(  this 
method,  and  it  is  the  purpose  of  this  paper  to  expose  an 
opportunity  to  use  this  idea  to  extend  Liapunov's  direct  method 
for  studying  stability.  Specifically,  wo  shall  be  interested  in 
extending  the  LaSalle  Invariance  Principle  [12]  to  the  context  of 
non a ut on omous  functional  differential  equations  (FDE) . 

We  ben  in  our  discussion  bv  giving  the  historical  development 
of  invariance  principles  in  Chapter  IT.  Then,  in  Chapter  111,  vo 
construct  a semi  flow  from  the  solutions  of  nonautonoinous  func- 
tional differential  equations.  From  this  semi flow  an  in- 
variance principle  can  be  obtained.  In  Chapters  IV  and  V we  weaken 
the  requirements  on  the  vector  field  necessary  for  obtaining  an 
invariance  propertv.  First,  in  Chapter  IV,  v;c  use  an  extension  of 
Ihc  concept  of  differential  equation  due  to  Kurxwe.il  [31].  The 
use  of  Kurxwe.il. ' s concept  of  a differential  equation  allows  us  to 
consider  the  stability  problem  for  a class  of  FDE  larger  than  that 
liscussed  in  Chapter  TU.  In  Chapter  V the  invariance  property  is 
obtained  directly  from  a continuous  dependence  result  without 
formally- constructing  a semiflow.  This  allows  us  to  consider 
equations  with  nonunique  solutions.  Also,  in  Chapter  V,  we  intro- 
duce the  integral- like  operator  equations  of  Neustadt  |43] . These 
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t ’ operator  equations  allow  for  yet  another  enlargement  of  the  class 

of  PDE  for  which  an  invariance  principle  can  be  obtained. 

Chapter  VI  contains  results  on  uniform  asymptotic  stability,  total 
stability,  and  eventual  stability.  All  of  these  results  are 
shown  to  be  intimately  connected  to  the  stability  properties  of 
the  limiting  equations  of  FOE.  Examples  illustrating  the  stability 
theorems  are  given.  Finally,  in  Chapter  VIT  we  indicate  the 
extension  oT  the  above  insults  to  FOE  with  infinite  delay  and  to 
neutral.  FDE.  Appendix  A and  M summarize  results  from  topology  and 
topological  dynamics,  respectively,  that  arc  of  importance  for  the 
body  of  the  text. 

Throughout,  wo  shall  be  concerned  with  working  directly  with 
the  differential  equations.  Dafermos  (11.5]  has  studied  the 
stability  problem  by  working  directly  with  the  solutions. 

A comment  on  notation  is  in  order.  We  abuse  two  symbols.  The 
symbol  w variously  stands  for  the  m-limit  set  in  dynamical  systems 
and  for  the  maximal  right  hand  endpoint  of  the  domain  of  definition 
of  a solution  of  a differential  equation.  The  symbol  variously 
stands  for  the  closure  operator  in  topology  and  for  the  operation  of 
taking  a one  sided  limit  from  below.  In  either  case  the  particular  j 

usage  of  the  symbol  in  the  text  eliminates  any  ambiguity. 
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governed  by  linear  funciional  dieeerential  equations 

by 
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Abstract:  Aspects  of  the  approxima t ion  and  optimal  control  of  systems  governed 

by  linear  retarded  nonau t onomous  functional,  differential  equations  (FDE)  are 
considered.  First,  certain  FDE  are  shown  to  be  equivalent  to  corresponding 
abstract  ordinary  differential  equations  (ODE).  Next,  it  is  demonstrated  that 
these  abstract  ODE  may  be  approximated  by  difference  equations  in  finite 
dimensional  spaces.  The  optimal  control  problem  for  systems  governed  by  FDE  is 
then  reduced  to  a sequence  of  mathematical  programming  problems.  Finally, 
numerical  results  for  two  examples  are  presented  and  discussed. 


Ah'?  »rao»  <4  Al’I’IiOXI  MA  C I f'M  AMI)  OP  I I MM,  COMTK''!,  OP  I.IMIIAI'  limilHTAKY  SYSTIllir, 
by  Douglas  (Mir  is  t ian  Peber,  Ph.l).,  lirown  Univerri  i v,  'Minn  1070. 

Our  concern  in  this  i tivcsl  i »*,  1 1 ion  in  with  the  - 1 1 »|  «i  i in. » I ion  and  opt.  inn] 
control  of  systems  p.overned  |,y  linear  rot  irded  functional  di  f I '•rent  in  I equa- 
tions ( I I'll ) . In  eh  apt  «r  T wo  or;  t ab  I i sh  t h"  existence,  uniqueness  and  coni  in- 
nous  dependence  of  solution';  of  ll)b.  Wo  furtfiee  demons  1 ra  t e t hi  t cor  lain  1011 
are  equivalent  to  cot  ■»•«  ■ «?p<  mv  I i nr  a'  ■ s 1 1 i<  I > u "I  i n n • v d i * I • > n t i . i I • ■ u i a • i • -n  ( < 'I  ) 
.'•'trli  0011  are  also  Inown  as  abs  t raw  1 evolution  e.  pin  t i •■•ir  (Alb).  'banter  II 
details  the  manner  in  which  * Ivse  All!’.  may  bo  approx ima*  ed  by  difference  equ  i- 
l ions  in  finite-dimensional  spaces. 

The  optimal  control  problem  for  systems  jyv/erned  by  rbb  is  then  reduced 
to  a sequence  of  mathematical  pr  open  mini np  problems  in  chapter  Ilf.  In  chapter 
IV  we  discuss  numerical  results  for  two  systems,  havinr,  applied  standard  tech- 
nique:! oT  numerical  analysis  to  compute  the  solutions  of  the  approx  i mn  I i np, 
problems.  The  first  of  these  systems  was  chosen  for  i1  simplicity,  so  that 
an  analytical  solution  would  be  readily  available;  the  second  system  is  assoc- 
iated with  a biochemical  process. 


